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ABSTRACT 

Predicting the colors of Luminous Red Galaxies (LRGs) in the Sloan Digital Sky Survey (SDSS) has 
been a long-standing problem. The g,r,i colors of LRGs are inconsistent with stellar population 
models over the redshift range 0.1 < z < 0.7. The g — r colors in the models are on average redder 
than the data (of the order 0.1 mag) while the r — i colors in the models are bluer (by 0.05 mag) 
towards low redshift. Beyond redshift 0.4, the predicted r — i color becomes instead too red, while the 
predicted g — r agrees with the data. We provide a solution to this problem, through a combination 
of new astrophysics and a fundamental change to the stellar population modeling. We find that 
the use of the empirical library of Pickles (1998), in place of theoretical libraries based on model 
atmosphere calculations, modifies the evolutionary population synthesis predicted colors exactly in 
the way suggested by the data, i.e., gives a redder r — i color, and a bluer g — r color, in the observed 
frame at z = 0.1. The reason is a lower flux in the empirical libraries, with respect to the theoretical 
ones, in the wavelength range 5500 - 6500 A. The discrepancy increases with decreasing effective 
temperature independently of gravity. This result has general implications for a variety of studies 
from globular clusters to high-redshift galaxies. The astrophysical part of our solution regards the 
composition of the stellar populations of these massive Luminous Red Galaxies. We find that on top 
of the previous effect one needs to consider a model in which ^ 3% of the stellar mass is in old 
metal-poor stars. Other solutions such as substantial blue Horizontal Branch at high metallicity or 
young stellar populations can be ruled out by the data. Our new model provides a better fit to the 
g — r and r — i colors of LRGs and gives new insight into the formation histories of these most massive 
galaxies. Our model will also improve the k- and evolutionary corrections for LRGs which are critical 
for fully exploiting present and future galaxy surveys. 
Subject headings: stars: HB galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Age-dating the stellar populations of galaxies pro- 
vides astronomers with a cosmic timescale which, being 
ruled by stellar evolution, is independent of cosmologi- 
cal models and allows the use of gala xies as cosmologi- 
cal probes (e.g. iJimenez fc Loebl '2002). The interpreta- 
tion of galaxy spectra in terms of stellar populations is 
also the only effective way of reconstructing their star- 
formation histories, which allows one to get clues on the 
still poorly-known process of galaxy formation. There 
is a long history of research in thi s area. We provide 
here a selection of recent works ( Bernardi et al.l l2006l: 
Jimenez et al. 2007; Kauffmann ct al. 2003; Nelan et al.l 
2005i : .Panter et al...2007; .Thomas et al. .,2005( 1 and refer 
the reader to Renzini (2006) for a comprehensive review 
of studies at both low and high redshift. 

Studies of galaxy evolution based on the Sloan Digi- 
tal Sky Survey (SDSS) data of Luminous (massive) Red 
Galaxies (LRGs) have highlighted a potentially funda- 
mental problem with the rest-frame optical of stellar 
population models. As first noticed bv l Eisenstein et aTl 
UpOl), the models are systematically too red in the 
g — r observed-frame in the redshift range 0.1 to 0.4 (see 
Figure 1, left-hand panel). This discrepancy could not be 
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cured by adopting more complex stellar population mod- 
els with various star formation histories, different stellar 
population model codes or empirical galaxy spectra. To 
be able to analyse the data, Eisenst ein et al.. (2001) ap- 
ph ed a shift of 0.08 mag to the models. 

(|Wake et al.l2006l W06) extended this study by includ- 
ing the r — i color as a further constraint and by adding 
the 2dF SDSS LRG and Quasar (2SLAQ) survey sam- 
ple Ic annon et al.l fe006l to increase the redshift range 
to z ~ 0.8, demonstrating further discrepancies between 
the models and data. While the models are too red in the 
observed frame g — r ai low redshift, the r — i color turns 
out to be too blue (see Figure 1, left-hand panel). This 
pattern changes with increasing redshift such that the 
predicted r — i color becomes too red, while the predicted 
g — r matches the data (all observed frame). Figure 1 in 
W06 demonstrates that the addition of star formation 
does not solve the problem. Indeed, the light of young, 
blue stars while curing the g — r syndrome at the low- 
est redshifts, worsens the discrepancy in the r — i color 
at low redshifts and in the g — r color at high redshifts, 
such that colors in the model become too blue. The ori- 
gin of this mismatch has remained a puzzle so far, and 
has serious consequences for the interpretation of current 
observations and the planning of future galaxy surveys. 
In this paper we provide a solution to this problem and 
present a working stellar population model. 

The paper is arranged as follows. In Section 2 we recall 
the features of the luminous red galaxy sample. In Sec- 
tion 3 we present the working solution and we comment 
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on other discarded options. A conclusive discussion is 
placed in Section 4. 

2. THE DATA SAMPLE 

W06 extract their sample o f LRGs from th e SDSS 
multi-color photometry fSDSS- l^k et al.|[2000l ). Spec- 
troscopic redshifts are taken from the SDSS in the red- 
shift range 0.15 < z < 0.37, and from the 2SLAQ survey 
((Cannon ct al. 2006) in the redshift range 0.45 < z < 0.8. 
At z > 0.4 we use the deeper multi-epoch data in stripe 
82. LRGs from both surveys were selected in a consistent 
way using the same rest-frame color cuts, under the as- 
sumption that they were old, passively evolving galaxies, 
in order to make sure that the same galaxy population 
is sampled at all redshifts (see W06 for details). 

Obviously, the sample definition depends on the color 
selection, which in turn depends on the model adopted 
to identify the LRGs. This is what we should call model- 
dependent data setting. One of the worse consequences 
of the discrepancy between data and models is that in 
order to minimise the effect of incorrectly /c-correcting 
the data a conspicuous part of the initial data base had 
to be discarded. This fraction was as substantial as 85% 
in W06. Hence, besides the general astrophysical goal of 
understanding the stellar populations of luminous, mas- 
sive galaxies, there is the more practical need to opti- 
mally exploit the precious data emerging from present 
and future surveys of LRGs to high redshift. 

To ensure self-consistency, the data used in this paper 
have been re-selected using the ()Marastonl 12001 M05) 
models. 

In the following we assume that galaxies start form- 
ing stars at a redshift of five, a s suggested by the lo- 
cal fossils (e.g. ,B crnarcl .ret al.ll200 6; Jimenez et al. 200'^ 
Kauffmann et al.ll2003t iNelan et al.i2005l : iThomas et all 
20051 ) ■ which fixes their present age to be about 12 Gyr^. 

3. THE WORKING MODEL 

The problem as explained in Section 1 is illustrated in 
the left-hand panel of Figure [Tl which is a remake of Fig- 
ure 1 by W06. The red line refers to a single-burst model 
(Simple Stellar Population, SSP) with solar metallicity 
and an age of 12 Gyr at redshift zero. Residuals are 
shown in Figure [2l The right-hand panels of Figures [1] 
and [2] show the solution to the problem using our new 
model. The latter is based on two modifications; the re- 
placement of theoretical stellar atmospheres by empirical 
ones, and the inclusion of a small (3% by mass), metal- 
poor subcomponent. Details of this model are described 
in the following sections. 

3.1. Inclusion of empirical stellar atmospheres 

The simultaneous mismatch in opposite directions of 
colors sampling close rest-frame wavelengths was sug- 
gestive to us of a shortcoming in the stellar spectra 
that are used in th e population synthesis models (e.g. 
iLeieune et al.lfl997l based on the Kurucz 1979 hbrary). 
To explore this path we consider empirical spectral li- 
braries as a substitute to the theoretical ones and com- 
pute stellar population models using such libraries as in- 
put to the M05 code. 

^ Wake et al. used 13 Gyr, but this small difference does not 
impact on the overall conclusions. 



Several high-to-intermediate resolution libraries of 
flux calibrated empirical stellar spectra, intended 
for use in evolutionary population synthesis , have 
been constructed o ver the last deca d e, e.g . iPicklesI 
(fl99l. STELIB (iLe Borgne eiTaH |2003), MILES 
dSanchez-Bla zQ ucz eTal.! I2006D and ELODIE.vS.l 
( Prugniel et al. r' l2007fl . For the aim of this work, 
which is to study the redshift evolution of spectra, 
the Pickles library is the only suitable one, due to its 
wide wavelength coverage (1150 - 25000 A). Models 
constructed from the other empirical libraries were also 
considered, but their wavelength coverage is too narrow 
(3900-6800 A for the EL0DIE.v3.1. and 3500-7500 A for 
MILES, respectively) to make them suitable for studying 
the evolution of spectra with redshift. 

The mapping of stellar atmospheric parameters in the 
Pickles library is very good, especially the lower main 
sequence and the tip of the red giant branch. The rather 
low spectral resolution of the Pickles library {R « 500) 
is not an issue here, since we are only interested in 
broadband colors. However, the library does not contain 
enough stars to cover all important evolutionary phases 
for non-solar metallicities. Therefore, this exercise fo- 
cuses on models with solar metallicity, the most relevant 
for massive galaxies such as LRGs. 

The effect of the input stellar library is shown in Fig- 
ure [3l which displays the spectral energy distributions of 
two 12 Gyr, solar metallicity SSPs models that only differ 
in their input spectral library. The red line shows the one 
based on the P98 empirical library, whereas the black one 
is based on the Lcicune ct al. ( 1997) theoretical library. 
A flux excess is evident in the theoretical spectrum be- 
tween ~ 5500 and ~ 6500 A. This discrepancy increases 
with decreasing effective temperature independently of 
gravity, and points to problems in the model atmosphere 
calculations at temperatures below ~ 6000 K. 

Also shown in the figure are the central wavelengths 
of the g, r, i SDSS filters at redshifts zero and 0.1. The 
lower flux in the empirical spectrum is sampled by the 
r- filter at redshift 0.1, which explains the redder r — i 
colors (by ^ 0.06 mag) of the empirically-based models 
with respect to the theoretical ones. For the same reason, 
the g — r at z = 0.1 is bluer by ~ 0.05 mag. The com- 
parison with SDSS data is shown by the middle panels 
of Figures [Hand [21 The inclusion of the empirical library 
rectifies the position of the models with respect to the 
data at the lowest sampled redshift. The predicted g — r 
color becomes bluer and the r — i color redder around 
z ^ 0.1. The trend is the opposite at higher redshifts; 
the observed-frame g — r color gets redder and the r — i 
color slightly bluer around z ~ 0.6. Hence the model 
has improved mainly at lower redshifts, and most signif- 
icantly in the r — i color. Further modifications to the 
model are still necessary to track correctly the evolution 
with redshift. However, the situation has become much 
clearer with the inclusion of empirical atmospheres; the 
model spectral energy distributions (SEDs) need to be- 
come systematically bluer shortward of the wavelength 
range sampled by r — z at z ^ 0.1. 

3.2. The metal-poor subcomponent 

As the middle panel of Figure [1] shows, the match be- 
tween models and data still requires a slightly bluer g — r 
colors at low redshift, but now both bluer g — r and 
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Fig. 1. — The g — r and r — i colors of LRGs as functions of redshift (points; data from W06). The median is given by the green Hne. 
Typical errors as function of redshift indicated by the error bars. Left panels. A solar-metallicity passively evolving, single-burst model 
with an age of 12 Gyr at redshift zero (red line). Middle panels. Same data as in the left-hand panel. The stellar population model uses the 
Pickles (1998) empirical spectral library instead of the theoretical one (see text). Right panels. Same data as in the left-hand panel with 
a composite model with 3% by mass of metal-poor stars. Both the metal-rich and the metal-poor component are 12 Gyr old at redshift 
zero. The metal-rich component uses the Pickles (1998) empirical spectral library. 
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Fig. 2. — Residuals of the plots shown in Figure [T] 



r — i colors at high redshift. The flatness of g — r be- 
yond z — 0.4 further suggests that the required bluen- 
ing component must be slowly evolving with look-back 
time. This disfavors the option of adding a component 
of residual star formation (see W06 and Section 3.3). A 
metal-poor old sub-component is the best candidate. 

Metal-poor stellar populations have blue turnoff 's and 
very often Blue Horizontal Branches (BHBs). A BHB 
from a metal-poor stellar population arises due to the rel- 
atively high effective temperature of the evolving stars, 
even without assuming large amount of mass-loss dur- 
ing the Red Giant Branch phase. Calculations show 
that a stellar population with metallicity [Z/H] ~ — 2.2 
has already developed a BHB at ages of 6 Gyr with 
a mass-loss efficiency consistent with calibrations from 
Milky Way globular clusters (Figure 11 in M05.) 

The difference in spectral energy distributions be- 
tween a solar-metallicity model and the same contam- 
inated by metal-poor stars (3% in mass) is such that 
the composite SED has more flux in the optical bands, 



while leaving unaltered the flux longward ~ 9000 A 
(|Maraston fc Thoma!j|2000f) . This helps the g - r color 
to become bluer without perturbing the r — i color at the 
lowest redshifts. At the highest redshifts it makes both 
colors bluer as required by the data. 

As mention above, components of metallicities other 
than solar cannot be based on the empirical spectral li- 
brary of Pickles. The metal-poor subcomponent is there- 
fore taken from the original M05 model based on theo- 
retical stellar atmospheres. Even if the inadequacies of 
the theoretical atmosphere in the optical wavelength re- 
gion at solar metallicities (see Figure |3]) propagates to 
the lowest metallicities, we can safely neglect this effect, 
as the metal-poor subcomponent contributes only 3% to 
the total mass of the composite stellar population. 

The final composite solution and the color residuals are 
shown in the right-hand panels of Figures [1] and [2l The 
redshift evolution of the SDSS data are now well matched 
by the model in both the g — r and r — i colors with a 
residual discrepancy of only 0.02 mag in r — j around 
z ~ 0.3. 

3.3. Discarded options 

The inclusion of a BHB in the metal-rich component 
is less attractive, as a very large fraction (50%) of the 
population would need to develop the BHB. 

We also tested the possibility of including recent star 
formation on top of the passively evolving SSP based on 
the Pickles stellar library. We included a small subcom- 
ponent with constant star formation rate since z — h 
summing up to a mass fraction of 3% at z = 0. We con- 
firm that the blueing of observed- frame g — r becomes too 
large around z = 0.6 with respect to z = 0.1. This effect 
could be avoided only by assuming that the level of resid- 
ual star formation in LRGs has been increasing steadily 
at least since redshift z — 0.7, which seems contrived. 

Finally, the LRGs in the SDSS are known to host a/Fe- 
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Fig. 3. — The effect of the input spectral library. The spectral energy distributions of 12-Gyr-old, solar-metallicity SSPs using theoretical 
(black) and empirical (red, from Pickles 1998) spectral library. Note the flux excess in the theoretical spectrum between ~ 5500 and 
~ 6500 A. Also shown are the central wavelengths of the g, r, i SDSS filters at redshift zero and 0.1. 



enhanced stellar populations ijEisenstein et al.l I2003D . 
We tested the effect of including such enhancements by 
calculating models using the a/Fe enhanced tracks of 
[Salasnich ct al. (2000). This yields a slight bluening of 
both g — r and r — i only at the lowest redshifts, while 
both colors at z > 0.4 remain unchanged. This is be- 
cause a/ Fe enhancement impa cts on the tracks only at 
old ages (jSalasnich et al.l [20001 ). Hence this option does 
not represent a viable solution. 

4. DISCUSSION 

In this paper we provide a working stellar population 
model for luminous red galaxies (LRG), which provides 
a substantial improvement to a long-st anding mismatch 
with the data from th e SDSS survey (jEisenstein et al.l 
[200lt I Wake eFall 120061) . 

The new model adopts empirical stellar spectra from 
the library of Pickles (1998) in place of the theoretical 
ones at solar metallicity in the evolutionary population 
synthesis. We found that an excess flux around 6000 A 
rest-frame in the theoretical spectra was responsible for 
making the synthetic r — i color too blue at redshift 0.1, 
and the g — r color too red. This finding impacts on a 
variety of studies which involve the B — V and V — R rest- 
frame and might be the cause for a mismatch between 
stellar population models and Milky Way Globular Clus- 
ter data in S — V" as highlighted in M05. 

The new model further includes the addition of a small 
(3% by mass), metal-poor subcomponent, being coeval to 
the old and metal-rich dominant population. The metal- 
poor old subcomponent can better match the observed 
color and color evolution because of its slow evolution 
with redshift in comparison to a young sub-population. 
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whose color evolution is much stronger. This model 
was constrained using data in the redshift range 0.1 to 
~ 0.8. An early version of this model was used in Cool 
et al. (2008) for studying the luminosity function up to 
a redshift 0.8. Though the model allowed for a better 
match towards the highest redshifts, the synthetic g — r 
was still too red at the lowest redshifts (cfr. Figure 13 in 
Cool et al. 2008), and the r — i was not well recovered. 
This early model did not include the empirical spectra. 
For this reason a larger metallicity for the dominant pop- 
ulation needed to be assumed, and as a consequence a 
higher fraction by mass of metal-poor stars. The use of 
the empirical spectra has solved the remaining problems. 

For the discussion of the possible origins of a metal- 
poor subcomponent in ma ssive galaxies we refer to 
iMaraston k, ThomasI (|2000l ). It needs to be assessed 
whether the requirement that 3% by mass of stars in 
LRGs are very metal-poor is just the reflection of the 
metallicity gradient known to be present in massive 
galaxies or is indeed caused by the presence of a very 
metal-poor subcomponent. The latter could originate 
from accretion of metal-poor dwarf satellites during the 
evolution of the galaxy. Future investigations on this 
question will be very valuable. 

The models are available at www.maraston.eu. 
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